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We demonstrate the possibility of turbulent flow of electrons in graphene in the hydrodynamic 
region, by calculating the corresponding turbulent probability density function. This is used to 
calculate the contribution of the turbulent flow to the conductivity within a quantum Boltzmann 
approach. The dependence of the conductivity on the system parameters arising from the turbulent 
flow is very different from that due to scattering. 
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Graphene continues to offer unexpected possibilities for theoretical studies which are novel and can be experimentally 
tested. In recent work [l], a regime has been identified for the electric transport in an ideal single layer undoped 
graphene sheet at zero applied gate voltage, where a hydrodynamic description using a quantum Boltzmann equation 
is relevant. The ratio of the shear viscosity r] to the entropy density s in this region was shown to be very small Q. 
This makes graphene behave like a nearly perfect fluid. It was also suggested that turbulent electronic flow in this 
regime might be possible, since a hydrodynamic system with low viscosity is expected to have a quantum turbulent 
flow [5|. A turbulent regime is characterized by a large value for the Reynold's number, which for graphene at a given 
temperature T is given by @ 
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where v is the Fermi velocity, Ut yp denotes a typical velocity and L is the characteristic length scale for the velocity 
gradients. Furthermore it was suggested in Q that for Re ~ 10 2 , turbulent flow might be possible in graphene. Thus 
determining the turbulent probability density function for graphene and using it to calculate its contribution to the 
conductance in graphene is of interest. These two issues are addressed in their Letter. 

Quantum turbulence is an approximate way to understand scaling behaviour in a quantum fluid in a nonequilibrium 
state. The framework for understanding the similarities and differences between classical fluid turbulence and quantum 
fluid turbulence is based on a nonlinear Schrodinger equation model Q or a filament model The theoretical 
approach we use may be viewed as a generalization of the nonlinear Schrodinger equation approach. It has been 
stressed by Zakharov et al Q that turbulence is a physical phenomenon which is much wider in scope than simply 
fluid turbulence. A quantum turbulent state is possible for a complex quantum system with many degrees of freedom. 
It represents a highly excited state of the quantum system where energy is pumped in one scale and dissipated 
due to nonlincaritics present at a different scale. Graphene is a complex system described by a nonlinear quantum 
Hamiltonian given by Eqn. (1). Furthermore, recent work h as established that this system has a hydrodynamic 
regime with a high Reynold's number. In view of this we would like to follow the well established methods 0-01 to 
determine the turbulent PDF for graphene and calculate the contribution made by this component to conductance. 

Turbulent flows are characterized by a probability distribution function (PDF) with scaling properties and an 
associated nonvanishing energy flux. Quantum turbulent flows have been studied using the methods of quantum field 
theory, as such systems have many degrees of freedom exhibiting stochastic behaviour Q. Systems where the wave 
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motion is a recognized feature and turbulent flows are generated due to nonlinear interactions are known as weak 
wave turbulent systems Q . 

In this Letter we calculate a possible turbulent flow probability distribution function (TPDF) for graphene in the 
hydrodynamic regime treating it as a weak wave system. The wave motion comes from the free Dirac equation. The 
nonlinearity in graphene is due to the electron-electron Coulomb interaction, which is described by a term quartic in 
the field operators We show that these combined types of interactions admit a solution with a turbulent flow. 
Using the TPDF we then calculate the conductance due to this flow. We find that the contribution to the conductance 
due to turbulent flow has a structure different from that found using scattering (lj. 

In Ref. [l| it was shown that graphene near the Fermi points is described by a Hamiltonian given by 
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where a+i, a~i are the Fourier mode operators, i,j = 1, 4 for graphene corresponding to two valleys and two spins 
and A = (+,—). We also define k = |kj, which is the modulus of the two momentum k. The explicit expression for 
T\ 1 \ 2 \ 3 \ i (ki, k 2 , k3, k 4 ) is given in Ref. [if. For our discussion, we need to note that under scaling, 
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where /i is a scalar. 

A turbulent flow represents a steady state of the system in which, unlike a thermal equilibrium state, there exists 
a transport of a charge current. Such flows are expected to have scaling properties. To calculate the PDF for such a 
flow, we follow the method described in Ref. modified to deal with quantum masslcss Dirac system. The basic 
idea is to calculate 
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between two "in-states" |8[, which are eigenstates of Ho- Next we look for scaling solutions which have the property 
that 
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dt 

and have a nonvanishing associated flux, to be discussed later. From Eq. (4) we find that 
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We would like to obtain the stationary solutions in terms of the energy instead of the momenta. We shall further 
assume that the solutions are isotropic, so that both n k ^\ and the energy e depend only on the magnitude k of the 
momentum vector k. Following the procedure of Zakharov Q we change variables from k to e(k), where the latter 
denotes the energy and replace T by 
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where we have used the relation d 2 ki — kidkidfli. From Eq. (6) it can be easily seen that the quantity U is invariant 
under scaling. 



From Eq. (5) and Eq. (6) and using the free particle density, we get 
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We now analyze Eq. (7). The 5-function can be used to carry out the €3 integral. What remains is a region D of the 
(ei, e 2 ) plane. In addition we have a constraint that e 3 = e% + e 2 — e > 0. Following [9j, we divide D into four sectors 
as follows 

Di = {(ei,e 2 ) G D \ e 1 < e, e 2 < e} 
D 2 = {(ei,e 2 ) G D \ ei > e, e 2 > e} 
D 3 = {(ei,e 2 ) G D \ ei < e, e 2 < e} 

Z? 4 = {(ei,e 3 ) GD I ei > e, e 2 < e} (8) 

It was shown in 0] that the regions D 2 , D 3 , D4 can be mapped into D\ using the Zakharov transformations. For 
example D 2 with variables (e' l7 e' 2 ) can be mapped to D\ with variables (ei, e 2 ) using 
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and so on. Using these transformations and the ansatz n(e) = Ce -1 , we get 
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where y = 3x — 3. The factor within the square brackets in Eq.(10) can be made to vanish by choosing y = or 
y = 1, irrespective of other factors in the equation. It is possible to show (El [l0| that only the solution corresponding 
to y = or x = 1 has an associated nonzero flux and hence represents a turbulent state. This is the solution that we 
shall use. 

Our solution for the TPDF is thus given by 
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where for the moment we have set Vp — 1. We now need to determine C. To do this, we observe that the total 
number of turbulent particles per unit area, that is the turbulent number density Nt is given by 
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where e max — vpk max — jIcbT, where T is the background temperature of the system before the turbulent flow 
starts, ks is the Boltzmann constant and 7 << 1 is a constant, which is required for the hydrodynamic regime to be 
relevant. 

We now proceed to calculate the conductance using the quantum Boltzmann equation. For this purpose we introduce 
a time dependent electric field E(t), which is assumed to be small. We write the time dependent density 
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and = t?) the turbulent time independent PDF and determine fx by using the equation 
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We use the ansatz 
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We find 
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where we have ignored the terms higher order in E. It is now a simple matter to calculate the conductance. Taking 
the electric field to be in the x direction, we have 

<rr(u) = (19) 

where {< Ji >) — fx and ut is the contribution of the turbulent flow to the conductivity. Substituting for f\ and 
introducing back the Fermi velocity Vp, we get 

e 2 N T v F f ~/k B TL \ ( 1 \ 

where w^fc = jksT, 7 << 1 and L is the typical sample size. It may be noted that the spin valley degeneracy factor 
does not appear in the contribution to the conductivity coming from the turbulence flow. It is useful to compare the 
conductivity or obtained above with the expression for conductivity due to scattering ct sc obtained in [![, which is 
given by 

_ e 2 Nk B Tln2 
h —ihjj + nksTa 2 ' 

where a is the electron-electron interaction strength and N denotes the spin valley degeneracy in graphene. Comparing 
the expressions of ctt obtained here with that of a sc obtained in [l[ we see that they have very different dependence on 
the physical parameters of the system. In particular, ctt(w) is independent of the spin- valley degeneracy in graphene. 
The temperature dependence of <jt and ct sc arc also quite different. 

To summarize, we have calculated the contribution of the turbulent flow to the conductance in graphene. The 
dependence of the conductivity on the system parameters is very different from that described in recent calculations 
[l[ using scattering. Note that the temperature in our formula is a natural cutoff describing the hydrodynamic region. 
It represents the background equilibrium state temperature of the graphene sample. The factor Nt present in the 
expression of the conductance is assumed to be of the same order of magnitude as the uniform electron density in 
the system. This corresponds to the physical picture of a finite macroscopic fraction of the electron fluid exhibiting 
turbulent behaviour. 

The turbulent flow discussed here is a novel consequence of nonlinear electron-electron interactions in graphene. The 
turbulent flow found has superfluid like properties in the sense that the dispersion relation is linear in the momentum 
and the system has a macroscopic turbulent flow density. The linear dispersion implies that dissipation through 
scattering with a defect is only possible when the fluid velocity is greater than vp- The absence of dissipation implies 
that the electron flow is frictionless in the turbulent regime. 

As emphasized in [l[, in order to experimentally access this region, a large graphene sample is required. With 
a suitable sample, the contribution of the conductivity due to the turbulent flow might be identified through its 
particular dependence on the system parameters. 



[1] L. Fritz, J. Schmalian, M. Muller and S. Sachdev, Phys. Rev. B 78, 085416 (2008). 

[2] M. Muller et al, Phys. Rev. Lett. 103, 025301 (2009). 

[3] E. P. Gross and L. P. Pitaevsky, Phys. Rev. 106, 161 (1957). 

[4] K. W. Schwarz, Phys. Rev. B 31, 5782 (1985). 



■5 



[5] V. E. Zakharov, V. S. L'vov and G. Falkovich, Kolmogorov Spectra of Turbulence I, Springer- Verlag, Berlin, 1992. 
[6] Y. V. Lvov, R. Binder and A. C. Newell, Physica D 121, 317 (1998). 
[7] Y. V. Lvov and A. C. Newell, Phys. Rev. Lett. 84, 1894 (2000). 
[8] M. Rakowski and S. Sen, Phys. Rev. E 53, 586 (1996). 

[9] V. E. Zakharov, Zh. Eksp. Teor. Fiz. 51 688 (1966) [Sov. Phys. JETP 24, 455 (1967)]. 
[10] D. Sanyal and S. Sen, Ann. of Phys. 321, 1327 (2006). 



